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Abstract 


Laminar- turbulent transition in high speed boundary layers is a complicated 
problem which is still poorly understood, partly because of experimental 
ambiguities caused by operating in noisy wind tunnels. The NASA Langley 
experience with quiet tunnel design has been used to design a new kind of 
short duration quiet flow tunnel which can be constructed less expensively. 
Fabrication techniques have been investigated, and inviscid, boundary layer, 
and stability computer codes have been adapted for use in the nozzle design. 
Construction of such a facility seems feasible, at a reasonable cost. Two 
facilities have been proposed: a large one, with a quiet flow region large 
enough to study the end of transition, and a smaller and less expensive one, 
capable of studying low Reynolds number issues such as receptivity. Funding 
for either facility remains to be obtained, although key facility elements have 
been obtained and are being integrated into the existing Purdue supersonic 
facilities. 
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1 Introduction 

The quiet-flow wind tunnel concept developed at NASA Langley over the past 
fifteen years is a major development in the experimental study of high-speed 
boundary layer transition (see, e.g., [3]). Unfortunately, current designs are 
beyond the reach of most laboratory and university budgets, a difficulty 
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which limits the amount of research progress possible. A quiet-flow Ludwieg 
tube design holds the promise of reducing the expense to a level where uni- 
versities and other laboratories could contribute. This report summarizes 
progress on the design, which the author believes is sufficient to show that 
the approach is feasible. A more detailed discussion of the motivation for the 
facility and the research to be conducted there is contained in [27], 

This final report enlarges on information presented in the semi-annual 
progress report [26], although that report was written after the bulk of 
the research funded under this grant was completed. A full summary of 
the progress will be reported in a more accessible form at the 1991 AIAA 
Aerospace Planes Meeting [27]. Progress was also reported on at the Novem- 
ber 1990 meeting of the Fluid Dynamics Division of the American Physical 
Society in Ithaca, New York. 

Since Ludwieg tubes have been around for many years, and NASA Lan- 
gley has already established the feasibility of creating quiet-flow wind tun- 
nels, the major question to be addressed was the cost of the proposed facility. 
Cost estimates were obtained for major system components, and new designs 
which allowed fabrication at lower cost were developed. A large fraction of 
the facility cost comes from the fabrication of the highly polished quiet-flow 
supersonic nozzle. Methods for the design of this nozzle were studied at 
length in an attempt to find an effective but less expensive design. Since 
the Mach number and Reynolds number of any tunnel would have to depend 
on the particular interests of the sponsor, and since a sponsor for tunnel 
fabrication has not yet been found, a specific nozzle design has yet to be de- 
termined. However, the improvement of the design tools is nearly complete, 
and many specific nozzles of various types have been studied, to varying de- 
grees of completeness. Less expensive methods for nozzle fabrication have 
also been investigated, and a test specimen for the fabrication te chn iques 
has been fabricated. Progress has been sufficient to show that a quality facil- 
ity can be fabricated for a reasonable cost. Instrumentation and fabrication 
techniques are being further investigated through modification of the small 
Purdue supersonic wind tunnel in order, to achieve quiet flow at low Reynolds 
numbers. 

The general design and most cost estimates are discussed in the first sec- 
tion. Computation methods used there but not specifically discussed were 
taken from Pope and Goin [24]. The test section configuration and shape 
involves many special considerations and a large fraction of the facility cost. 
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Design methods for the test section are discussed in the second section. Fab- 
rication methods have a major impact on the test section cost and are dis- 
cussed in the third section. The fourth section discusses improvements which 
are being made to the Purdue Aeronautics blowdown supersonic wind tun- 
nel, for the purpose of low Reynolds number work and instrumentation tests. 
The report concludes with a summary. Appendices contain equipment cost 
estimates for the large and small facilities. 

2 Design Overview 

2.1 General Description 

The Ludwieg tube wind tunnel is a long pressurized tube with a supersonic 
nozzle on the end (see Figures 1, 2, and 3). When the quick-opening 
valve opens, fluid flows from the tube through the contraction, throat, and 
supersonic expansion, through the test section, past the second throat, and 
into a vacuum tank. After some startup time which depends on model size, 
valve opening time, and test section configuration, the test section flow is 
essentially steady. The rapid expansion of gets from the tube sends an expan- 
sion wave upstream into the tube. This wave reflects off the fan 1 end of the 
tube; on its return to the test section, the useful test time has ended. At this 
point, the tube and test section are still pressurized. Fluid continues to flow 
into the vacuum tank until atmospheric pressure is reached downstream, and 
then blows out the flapper valve into the atmosphere until the tube has de- 
pressurized. The vacuum tank allows runs to be made at low total pressures, 
and assists in starting higher pressure runs. Some of the many published 
studies of various forms of the Ludwieg tube can be found in [11], [17], [18], 
[16], [25], [29], and [31]. 

The tube will be used for quiet-flow study of boundary layer transition. 
Here, it is desirable to reach transition Reynolds numbers on a flat plate, 
at all Mach numbers for which heating is unnecessary. The thickness of 
the laminar boundary layer on a flat plate was estimated using compressible 
boundary layer similarity theory. The equations were rederived following 
White ([30]). A computer program was written to solve the boundary layer 
and isentropic expansion equations and to generate a table of test section 
Mach number, stagnation pressure and temperature, model length required 
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Figure 1: Sketch of Proposed Quiet-Flow Ludwieg Tube 
High Reynolds Number Version 
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Figure 2: Sketch of Proposed Quiet-Flow Ludwieg Tube 
Low Reynolds Number Version 
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to reach transition Reynolds number, and so on. This table was used to create 
various contour plots to optimize the design choices. It became obvious that 
the best way to design the tunnel was to make the largest possible test 
section, so that the model and tunnel wall boundary layers are as thick as 
possible. 

Maximum allowable tunnel pressure should be sufficient to reach flight 
transition Reynolds numbers at the highest Mach number of interest, with a 
plate that will fit in the quiet-flow part of the test section. Since quiet-flow 
nozzle design is non-trivial, only estimates can be made at this stage. If a 
quiet-flow length of about half a meter can be obtained, then a stagnation 
pressure of 150 psi is sufficient to reach the transition Reynolds number 1 
of 2.1 x 10 7 at Mach 4. Higher Match number work would require a larger 
quiet-flow test section or a higher pressure. The initial test section is to be 
15 inches wide, so that it can be machined in the Purdue Central Machine 
Shop numerically controlled mill. A reasonable test section height of about 
6 inches makes for a mass flow rate of about 10 kg/sec. 

The useful tube run time also depends on the length of time needed to 
establish the flow. Work by Johnson et al. [13] showed that as expected this 
starting time is severed times the time needed for a particle to cross the test 
section (i.e., roughly 10-20 milliseconds). Although workers in the AEDC 
tube struggled with a much longer startup time, this was due to the other 
issues involved with their transonic test section and its slotted wall (see [29]). 
Several other studies of the startup time have since been conducted (see [6], 
[14], [15], [21], [20], [33]). These articles also agree as to the general criterion 
for the length of the starting time (which does differ considerably if the nozzle 
is a vented transonic design, as opposed to the supersonic designs considered 
here). It appears from the work of these authors that a starting shock may or 
may not appear in the test section during the starting process, depending on 
geometry. If such a shock is present, it may or may not damage any hot wires 
which might be used for flow measurement. A determination of the ability 
of standard hot wires to withstand such shocks is currently being carried out 
at Purdue as part of another project. 

Extrapolated from flight data presented in [8, Figure 8]. 
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2,2 Sizing of Tube 

The length of the driver tube governs the useful flow duration; the longer 
the tube, the longer the useful flow 3 . The Aerospace Sciences Laboratory at 
Purdue University is a 250 foot long aircraft hangar which has been converted 
to a lab. If 50 feet is allowed for the test section, valve, diffuser, and so on, 
then there is room for a 200 foot long tube. Such a tube would have a runtime 
of about 350 milliseconds. In order to make such a long tube useful, it 
must have sufficient diameter. A boundary layer grows behind the expansion 
wave propagating into the tube (see [23] for a fairly recent discussion). The 
displacement thickness effect of this layer causes a pressure variation in the 
test section which must be small for good flow. An experimental correlation 
for this variation is given by Russell et al. [25, Equation 1]. This variation 
depends on the mass flow rate out of the tube, which in turn depends on the 
contraction ratio and the test section size and Mach number. It would be 
desirable to make the tube large enough to allow for the future use of larger 
test sections, so that thicker boundary layers can eventually be studied. 

Since the cost is fairly insensitive to tube diameter, and since pressure 
uniformity is important, fairly large diameters are currently planned. Plans 
for the large facility call for a 48 inch diameter tube, this being the largest 
standard carbon steel pipe size. If the first stage of the contraction is ma- 
chined axisymmetric on a lathe, the contraction cost is not excessive even for 
this large diameter. For a 6 x 25 inch wide Mach 4 nozzle at total pressure 
120 psi, this diameter gives a pressure drop of 5 x 10~ 4 %, allowing ample 
room for larger mass flow test sections. A somewhat smaller, higher pressure 
tube may eventually be used, to match the tube size to the outlet plumbing 
size, and to maximize test section Reynolds number; for a 24 inch tube and 
the same test section, the pressure drop is 0.03%. The smaller facility plans 
call for a tube made from standard 12 inch steel pipe, some of which we have 
on hand from the vacuum system for our small supersonic wind tunnel. If 
the tube is 75 feet long, the running time will be about 0.15 seconds, more 
than sufficient to study instability waves. For this tube and a 3 by 6 inch 
Mach 5 nozzle at 120 psi total pressure, the pressure drop is 7 x 10“ 4 %, also 
allowing room for larger test sections. 

The facility working pressures are being designed to allow the facility to 
reach the maximum Mach number possible without adding the complexity of 

2 Other methods of extending the run time [1] do not seem cost-effective. 


8 


driver gas heating, and to maximize the thickness of the boundary layers in 
the test section. A driver section maximum pressure of 150 psi seems to be 
a good compromise between future flexibility and current cost, since the use 
of higher pressures merely makes for thinner test section boundary layers. 
The precise working pressure will depend on the test section sizes, which are 
currently still under study. 

Plans call for operating the facility primarily at room temperature, to 
save costs. However, it turns out that a standard carbon steel pressure vessel 
can withstand some 200°C of heating without any special treatment or extra 
cost; this bonus gives the facility the potential to reach low hypersonic Mach 
numbers in the future, when a tube heating apparatus could be added 3 . De- 
pending on the pressure and the liquefaction computation, a Mach number of 
6 to 7 could be reached, which is more than sufficient to study the hypersonic 
second mode instability waves (see [24, Figure 1:39]). 

2.3 Contraction and Test Section 

Plans are to make the contraction for the large facility in two parts with a 
flanged joint. This will allow for varying test section size without replacement 
of the entire contraction, and the joint should not give trouble when placed 
reasonably far upstream. The contraction is to be cast from carbon steel 
about 1 inch thick using a special one-time casting technique (stryofoam 
mold), and then machined axisymmetric to a specified contour using a tracer 
lathe. A price estimate of $16,000 was obtained from Frankton Machine and 
Tool, Inc., of Indiana; a slightly higher quotation was obtained from another 
firm. This cost is insensitive to contraction thickness and also insensitive to 
diameter (reduction of diameter by factor 2 reduces cost by about 30%). An 
extra $9,000 has been budgeted for the extra costs involved in building the 
contraction to ASME code; extra costs will also be involved if the second part 
of the contraction involves a transition from axisymmetric to 2D. These costs 
will depend on the diameter at which the transition is made. The contraction 
for the smaller facility would be hogged out of solid al uminum . 

The detailed design of the test section is discussed in a later section. 
A general issue involves the choice of a 2D or axisymmetric test section. 

3 A 600 amp 250 volt DC motor generator set and two 30 volt 200 amp DC power 
supplies are available, to be used to heat the tube by using the tube itself as a luge 
resistor. 
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Current plans are for a 2D test section, which is easier to machine and polish, 
and which allows easy optical access. Nozzle housings can be designed which 
allow for interchangeable nozzle blocks for different Mach numbers. However, 
axisymmetric test sections do not have problems with side-wall boundary 
layer contamination or corner vortices. It might also be possible to make 
a mandrel for am axisymmetric test section on a diamond turning lathe and 
obtain a high-quality surface without polishing. Such am axisymmetric nozzle 
would have to be machined with tighter tolerances on the surface contour, 
to avoid difficulties with focusing of weak shocks on the centerline. 

2.4 Valve Location and Type 

For quiet test section flow, the valve must be located downstream of the test 
section. Otherwise, disturbances generated by the open valve will disturb the 
flow. This means that a large diameter valve must be used, for the tunnel 
flow area is smallest at the first throat, and much larger at the second throat 
and downstream. For a test section of sufficient size the downstream flow 
area corresponds to a diameter larger than 12 inches. The valve must open 
in a time the order of 10 milliseconds so that it does not significantly reduce 
the runtime. Mechanical valves of this type seem to be very expensive. Thus 
it seems preferable to use a burst diaphragm for this tunnel, just as is done 
in shock tubes. Current plans are to use a pair of burst diaphragms. The 
tunnel is pumped up to half pressure, and air is bled into the area between 
the diaphragms. Then the tunnel is fully pumped up. After air is bled into 
the region between the diaphragms, the second and then first will burst at 
a time controlled by the bleed time. This allows for more precise control of 
the tunnel total pressure and thus the test section Reynolds number. The 
burst diaphragm design could be adapted from those used on the Caltech 
shock tubes (drawings have been obtained) or from those used on double 
diaphragm systems at CALSPAN. A simple design used at UT Arlington 
[32] might also be used. 

If a mechanical valve capable of rapid closure were obtained, successive 
runs could be made without complete depressurization of the tube, saving 
on pumping costs and time. However, a valve with appropriate specifications 
(about 1 square foot of open area, opening time about 10-20 msec) seems 
unavailable without custom engineering at a prohibitive cost (upwards of 
$50,000). This option seems best reserved as a future possible upgrade to 
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the system, should it seem desirable; it would be easy to bolt in a new valve 
if one was fabricated or obtained. 

2.5 Second Throat, Diffuser, and Vacuum System 

A variable second throat allows for better pressure recovery in a supersonic 
wind tunnel, but for a Ludwieg tube the flow time is fixed by the tube length 
and not the pressure recovery. The cost of a larger vacuum t ank is much 
smaller than the cost of a variable second throat. Thus, we anticipate using 
a fixed second throat, with its geometry linked to the test section geometry. 

A conical diffuser is a relatively inexpensive way of getting pressure recov- 
ery downstream. Since test sections of varying size are envisioned, it seems 
best to make the diffuser in two parts, so that the first only can be varied for 
smaller test sections. This also allows for the use of one valve of fixed size 
for a range of test sections. 

The size of the downstream exit plumbing is a limitation on the size of the 
test section. This downstream plumbing can be made in Purdue University 
shops at relatively low cost, since it does not have to be ASME code stamped 
(not pressurized in normal operation). Large pipe flanges, in particular, have 
to be custom made; the largest pipe flange which the Central Shop is capable 
of making is for pipe of about 30 inch diameter. Thus, plans call for 30 inch 
diameter piping. 

The vacuum tank size is controlled by the run time and by the mass flow 
rate. Since the run time is short a small vacuum t ank can be used; this is a 
very large cost savings. A vacuum tank of 500 cubic feet has been procured 
for use in the upgraded Purdue Supersonic Wind Tunnel and this tank is 
fitted with a 30 inch diameter welded cap so that it can be readily hooked 
up to the outlet piping for the proposed large Ludwieg tube. This size is 
sufficient to run the large tube at about Mach 4 for the full 350 msec as long 
as the flow total pressure is at least about 5 psia (assuming the vacuum pump 
can bring the tank down to 0.01 psia). At lower Mach numbers the run-time 
would be reduced if the total pressures were this low. The tank also contains 
a 12 inch opening which is to be connected to the existing supersonic wind 
tunnel, and which could be connected to the small Ludwieg tube. This tank 
was procured using funding from this grant and the School of Aeronautics 
and Astronautics at a total cost of about $14,000, including footings. 

A Stokes Model 212-H 150 CFM vacuum pump has been retrieved from 
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storage and is being restored to running condition. This pump can pump the 
500 cubic foot tank down to 0.01 psia in about 35 minutes, and is capable of 
reaching an ultimate vacuum of 10 microns of mercury, which is much less 
than will be required. 

The flapper valve is necessary to exhaust the flow when the tube is de- 
pressurizing after the test is over and the vacuum tank is full, or, for higher 
pressures, for direct flow from the tube. Such a valve has been fabricated as 
a hinged cover to the end of a pipe tee. 

2.6 Compressor and Filter System 

Since the system will be constructed from scratch to be a quiet flow tunnel, 
the whole system can be maintained like a clean room. This allows the air 
to be filtered during the slow pump-up phase rather than during the rapid 
air flow testing phase, making the filters much cheaper. Appropriate filters 
can then be obtained for a few thousand dollars. 

A Van Air Model 350-HL twin-tower heatless dryer is present in the ex- 
isting system. A precision dewpoint sensor capable of measuring dewpoint to 
— 120°C' has also been obtained and installed. This sensor includes electron- 
ics which should make it easy to control the switching of the dryer towers in 
response to the gas humidity. 

The Aeronautics supersonic facilities include an Ingersoll-Rand PA50 
215SCFM 120 psig compressor which would also be used for the Ludwieg 
tube. This compressor is capable of pumping the the 4 foot diameter 200 
foot long driver tube (2000 cubic feet) up to 135 psia in about 200 minutes, 
allowing for several runs in a normal working day. A second stage compressor 
would have to be added for higher pressure work, at higher Reynolds num- 
bers. A matched oil-free second stage compressor capable of reaching 335 
psig at the same flow rate was priced from Corken International at about 
$14,000; a lower pressure unit would be less expensive. This item has not 
been included in the budget, for it currently does not seem essential to the 
initial plans. 


2.7 Safety Issues 

Purdue University safety office personnel (Mr. Mike Kopas and others) have 
been included in planning from the early stages. All Purdue pressure vessels 
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must be in accordance with the rules established by the Indiana State Boiler 
and Pressure Vessel Rules Board. Contact has also been established with 
this Board, through the secretary Mr. Bud Meiring. One of the Board 
members is a Purdue faculty member, Prof. Jim Hamilton, which facilitates 
communication. An ad hoc committee of the School of Aeronautics and 
Astronautics has been formed to consider facility safety issues and has been 
kept aware of the facility safety issues. 

Since the tube is to be an approved pressure vessel, the Safety Office 
does not see a problem with the tube being in the same room as the ma- 
chinists, students, and other building occupants. The price estimate for the 
tube includes fabrication and installation by an Indiana State approved pres- 
sure vessel manufacturer in accordance with the ASME Pressure Vessel Code 
(section VIII, division I). It appears that the contraction will also have to 
be fabricated and code stamped by an approved manufacturer; this extra 
cost has been estimated in the budget also. The test section, initial diffuser, 
and valve sections are the only other sections which are pressurized for long 
periods of time to high pressures. Since it would be awkward to have these 
pass through the hands of a third party manufacturer, the plan is to de- 
sign these in accordance with the Code and have the designs checked by a 
professional engineer familiar with the relevant codes and standards. These 
sections would then be fabricated by the Purdue Central Machine Shop and 
hydraulically tested. It should be noted that ASME Code approved welders 
are available on campus. The sections would then be approved and operated 
as Indiana State Special pressure vessels. This procedure is not unusual for a 
university facility and will probably be necessary since vendors interested in 
fabricating the precision machined test section have not been found. Plans 
call for designing the test section so that the pressure containment vessel 
is independent of the supersonic nozzles, as in the Soviet supersonic wind 
tunnels at Novosibirsk. 

The tube components downstream of the valve are not normally pressur- 
ized (although they will sustain some pressure during operation) and thus 
need to conform to the ASME Piping Code rather than the Pressure Vessel 
Code. No special fabrication stamps are then required. These sections are 
to be fabricated from standard steel piping by the Purdue Central Machine 
shop and other Purdue machinists and welders. Copies of the relevant por- 
tions of the ASME Pressure Vessel and Piping Codes have been obtained for 
assistance in the design. 


13 



The driver tube supports and the test section working platform will be 
designed and installed by Purdue internal departments to usual standards ac- 
ceptable to the Safety Office. Estimates for these components were obtained 
and are presented in the budget. 

Thus, the only safety issue which appears to impact the cost of the facility 
is the requirement for the tube to accord with the pressure vessel codes. The 
long driver tube will be contracted out to a pressure vessel manufacturer, 
following the usual procedure. The smaller test section area components will 
probably be designed and fabricated in house, and then pressure tested. Con- 
siderable consultation haw revealed no other major safety related expenses. 

2.8 Other Issues 

The Ludwieg tube will make a considerable amount of noise while the flow is 
dumping to atmospheric. Fortunately, the Aerospace Sciences Lab is located 
in an area where the creation of loud noises is acceptable. The lab is at 
the Purdue University Airport, within 100 yards of the end of the principal 
runway. Regularly scheduled propeller planes and occasional jets land very 
close, so the added noise will not be all that noticeable. There is no private 
land nearby, and the nearest student housing is about a quarter mile away. 
The noise from the tube blowdown can be roughly estimated using results 
obtained by Starr ([29]). For a somewhat different configuration, he gives 
data showing that the noise would be about lOOdB at 200 feet, without any 
muffling, which is the limit given for residential areas. Since the nearest res- 
idential areas are much further away, a minimum amount of muffling should 
make the flow acceptable. It may be desirable to add more muffling later 
if discomfort to operators and building occupants is large; this seems best 
determined after installation. 

Although the Aerospace Sciences Lab has the 250 foot length needed to 
contain the large tube, there is not enough floor space for it. Fortunately, 
it is 18 feet from the floor to the building rafters, so there is plenty of room 
for supporting the tube above the floor. Plans call for 12 foot columns 
spaced every 20 feet to support the tube. The building columns are also 
spaced 20 feet apart and will be tied into the support columns to stabilize 
the tube laterally. A 20 by 50 foot mezzanine will be built using standard 
steel mesh platforms to provide a floor at a convenient height for the area 
around the tube test section. This mezzanine will also include posts necessary 
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to support the test section. Furthermore, a one ton trolley and hoist will be 
installed to allow easy installation and modification of the test section area 
tunnel elements and models. These building modifications will be carried 
out by Purdue internal departments and the costs presented in the budget 
represent formal estimates prepared by them. The small facility would need 
to be raised only about 7 feet above floor level, and a small working platform 
could easily be constructed. 

Finally, the reaction loads from the tunnel flow need to be accounted for. 
The maximum force can be conservatively estimated from the sonic velocity 
at the nozzle throat times the maximum mass flow rate there. For the initial 
test section for the large facility, this comes out to be about 1000 pounds of 
force. Thus, reaction load bracing can be limited to guy wires fastened to 
the floor. If larger test sections are eventually constructed more elaborate 
reaction load bracing may have to be constructed. 

3 Quiet-Flow Supersonic Nozzle Design Meth- 
ods 

The author had the privilege of spending eight weeks of the summer of 1990 
studying quiet tunnel design at NASA Langley, which has the only high 
Reynolds number quiet facility in the world. The main purpose of this trip 
was to learn the design methods developed by the lead engineer, Ivan Beck- 
with, and his coworkers, which include nozzle designer Frank Chen. Although 
the author suggested and implemented several modifications to the existing 
methods the framework of this discussion is that of the Langley quiet tunnel 
design methods (see Figure 4). 

Inherent to the idea of a supersonic nozzle is the design of the walls, 
following inviscid supersonic flow theory, in order to produce a shock free 
isentropic expansion to a uniform parallel flow, which is the usual test section 
requirement. This nozzle design is complicated by the requirement for a 
suction slot upstream of the nozzle throat to suck off the contraction wall 
boundary layer. This inviscid part will be discussed in the first subsection. 
Since the test section necessarily has viscous boundary layers, these are often 
computed to allow for a correction to the nozzle wall shape. This correction 
is discussed in the second subsection. Finally, for a quiet flow test section, it 
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Figure 4: Quiet Flow Nozzle Design Schematic 


16 


is important that the nozzle wall boundary layers be kept laminar as long as 
possible. The stability and transition of the boundary layers can be estimated 
using roughness estimates and e N theory discussed in the third subsection. 

The design method current when the author arrived at Langley in June 
1990 involved the use of three separate computer codes, ran on the old Lan- 
gley NOS Cyber 205 machines. These machines have primitive operating 
systems, and the use of the codes involved a considerable amount of data 
file editing using primitive editors. The complete analysis of a single noz- 
zle shape involved a great deal of operator intervention, and a considerable 
amount of waiting - a matter of weeks was involved. 

The author has automated the use of the three large computer codes 
required for the design process. This streamlining was achieved by adapt- 
ing the design codes into modern FORTRAN-77 (from FORTRAN IV) and 
porting them to modern machines. The codes are run separately, as before, 
but special output files are written from each code in a form suitable to be 
read by a simple interface program, which produces files that can be used 
as the input file to the next program in the chain. This scheme will allow 
the complete series of codes to be run automatically on a specified nozzle 
shape through use of a command or batch file. Thus, several nozzle shapes 
can be investigated in the course of a single night’s computer ran, instead 
of several weeks of computation and editing. The scheme is not complete, 
due to problems encountered in upgrading from the old stability program 
GORTLER to the new version E**MALIK. However, it has been successful 
so far, and promises to make the design process faster and simpler. 

3.1 Inviscid Compressible Design 

Inviscid supersonic nozzle design is not yet a standardized procedure. A 
good recent tutorial is contained in the textbook by Zucrow and Hoffman 
[34, Sections 15-5 and 16-4]. The author has had the benefit of several long 
discussions with Professor Hoffman, also at Purdue, who specializes in nozzle 
design, although of the rocket variety. The supersonic flow is hyperbolic, so 
that downstream conditions are set by conditions upstream, and the flow is 
computed using the method of characteristics. However, the boundary con- 
ditions to be used are a matter of design judgement. Computations normally 
begin in the nozzle throat, with the best computations using a transonic per- 
turbation scheme to compute the flow near the throat, assuming the flow is 
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nearly parallel there. Thus, the upstream subsonic flow must deliver a nearly 
parallel flow in the throat. These transonic perturbation schemes are only 
valid for Mach numbers very near 1, and require as input some information 
regarding the shape of the nozzle near the throat, usually in the form of the 
throat radius of curvature (the higher order terms usually being neglected). 

For earlier NASA Langley designs these transonic perturbation approxi- 
mations for the throat have been extended upstream to find the inner contour 
of the boundary layer bleed lip. However, the only requirement for this inner 
bleed lip is that it deliver parallel flow to the nozzle throat. The designer 
thus has a range of choices for this upstream contour which can be used to 
simplify the mechanical and structural design. The outer side of the bleed 
lip has in the past been designed to simple curves, requiring only that the 
flow not be turned too dramatically, and that the bleed slot contain a sonic 
region to reduce the amount of noise which can propagate into the test sec- 
tion. This seems reasonable, give the limited number of requirements on the 
slot geometry. 

Besides the inputs required for this transonic calculation, the nozzle de- 
sign also requires some further inputs. This further input can take at least 
two forms: First, the designer can specify the distribution of Mach number 
along the nozzle centerline for some distance (until the downstream parallel 
flow requirement takes over), or second, the designer can specify the initial 
shape of the nozzle wall (again, until the downstream parallel flow require- 
ment takes over). Thus, an inviscid nozzle computation requires the designer 
to specify conditions in the transonic portion of the nozzle and in the initial 
supersonic region. These requirements are in addition to the specification of 
parallel exit flow and exit flow Mach number. 

The design method used for existing Langley quiet nozzles 4 involved an 
inviscid flow code adapted by Frank Chen from the Nelms minimum length 
nozzle code [22]. The Nelms code was originally written to design supersonic 
rocket nozzles for minimum length. This kind of design involves the use of 
a sharp corner expansion in order to produce the minimum length. Chen’s 
code used Hopkins and Hill’s perturbation technique [12] for computing the 
transonic flow in the throat, at least for the axisymmetric nozzles. The sharp 
nozzle corner, inappropriate for a wind tunnel nozzle, was avoided by using 

4 See, e.g., [7]. Frank Chen has not written any detailed description of the nozzle design 
procedure. Descriptions of existing procedure are based on discussions with Frank Chen. 
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for the nozzle contour one of the inner inviscid streamlines; this technique 
is discussed in Nelms’ paper. The particular streamline chosen was unclear. 
The technique was modified in the late 1980’s to add a region of radial flow 
between the initial expansion region and the region where the wall is shaped 
to turn the initial characteristics so as to produce a uniform exit flow. This 
radial flow region allows the boundary layer to grow without any concave 
curvature and was thought to reduce the Gortler instability problem. 

Since Chen’s source code was unavailable, not documented, and ran only 
on the old NOS machines, an alternative was sought. A search of avail- 
able codes unearthed the Sivells design code [28]. This code was specifically 
designed for production of wind tunnel nozzles, and incorporated various 
special considerations to improve the uniformity of the flow. The program 
was well-documented, reasonably well written, and source code was available. 
The program computes both 2D and axisymmetric nozzles, and its author 
was recommended by Ivan Beckwith as a person who did careful work and 
produced good quality nozzle designs. Furthermore, this program allowed 
for the use of a variable region of radial flow, just as had been incorporated 
into the Chen code (see [7]). 

The Sivells code was acquired with the aid of Charles Johnson of NASA 
Langley, and adapted to run in FORTRAN-77 on an IBM AT clone. The 
code allows the use of a simple turbulent boundary layer computation scheme, 
which is not used. It incorporates the Hopkins and Hill transonic flow scheme, 
which has also been adapted for use in 2D nozzles in a carefully documented 
way. Upstream nozzle conditions are specified through specification of the 
Mach number distribution on the nozzle centerline. This distribution is kept 
continuous to keep continuous second derivatives in the nozzle wall shape, a 
condition which may be required in order to achieve smooth flow. The free 
parameters which control the nozzle shape are easily set, and the program 
runs on an 8MHz IBM AT clone in a few minutes, allowing a large number of 
nozzle shapes to be easily investigated. A subroutine was added to the code 
to print a special output file of the exact form needed for the boundary layer 
computation code, so that rapid computation of all the nozzle parameters 
can be achieved in a batch file without operator intervention. 
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3.2 Laminar Viscous Boundary Layer Computation 

The existing Langley design method used a FORTRAN-IV version of the 
code written by Harris and Blanchard [10] to compute the viscous boundary 
layer. This code had been heavily modified to change the output form and to 
produce wall radii of curvature information for the stability computations. 
Since this made debugging and testing the code difficult, it was decided 
to get a current FORTRAN- 77 version of the code direct from Harris and 
Blanchard (through Venkit Iyer at Langley), and write a separate program 
to take the standard output form and specialize it. This was done. The 
code runs on the IBM AT clone in about 30 minutes, with a reasonable grid, 
and runs much faster on a bigger machine. It has been tested on one of 
Harris’ standard test cases (in fact, the author has adapted the code to the 
department mainframe and given it to Purdue undergraduates, who use it 
for a class project where they compare turbulent boundary layer results to 
experiment). A new subroutine was also written for this code to produce 
output in the right format for ready conversion to the input format for the 
compressible stability code E**MALIK. 

3.3 Computation of Laminar Boundary Layer Insta- 
bility 

The test section design requires an estimate of the position of boundary layer 
transition on the tunnel walls. Sound is radiated downstream along Mach 
lines from the initial location of transition and from the turbulent bound- 
ary layer downstream. This sound contaminates the flow at all positions 
downstream of the Mach line from the location of transition. The length of 
the quiet-flow test core in the test section is determined by the streamwise 
distance between the beginning of uniform flow and the end of quiet flow. 
Design of a test section for maximum length of quiet flow thus depends on 
moving the transition as far downstream as possible. 

The existing version of the Langley nozzle design methods involved the 
use of the GORTLER instability code written by Malik, which computed 
Gortler instability on the nozzle walls. This was used since it was discovered 
that Gortler instability was primary for the designs usually used. The code 
computes the maximum growth of Gortler disturbances, using an e N tech- 
nique. However, this code involved a great deal of operator interaction, since 
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only one Gortler wavelength could be tested at one time. It was decided to use 
instead the new version of Malik’s instability codes, E**MALIK, which com- 
putes Tollmien-Schlicting type instability as well as Gortler instability, runs 
more automatically, and is written in FORTRAN- 77 instead of FORTRAN- 
IV. Dr. Malik graciously supplied the source code for this program, which, 
however, can still only run one Gortler wavelength at a time. Comments and 
write statements were added to the code to make it more user-friendly, and 
the i/o file structure was also modified, for the same reason. The code was 
also adapted to run multiple Gortler wavelengths in a single use. This code 
was run on the Langley Convex machines, and successfully reproduced test 
case 6 from a paper published by Malik [19]. These same results could also 
be reproduced using results transferred from the boundary layer code which 
is used for the nozzle designs. However, the author had difficulty getting the 
code to work on sample nozzle test cases. During the last week of his sum- 
mer 1990 stay at Langley he was helped by another user, who also informed 
him that there was an updated version, which was free of these bugs. This 
updated version has been obtained and is now r unning on a Sun Sparcstation 
2 at the Purdue Aerospace Sciences lab, but the author has not yet applied 
it successfully to the nozzle problem. 


4 Quiet-Flow Supersonic Nozzle Fabrication 
Methods 

The crucial issue in quiet-flow supersonic nozzle design is the delay of tran- 
sition relative to the initial location of the uniform flow region, in order to 
maximize the quiet-flow test core. Besides contouring the nozzle to tailor the 
pressure gradient and curvature in order to reduce the growth of instability 
waves, it is also necessary to smooth the tunnel walls so that transition is 
not tripped by small roughness elements, which can be large compaxed to 
the thin accelerating supersonic boundary layer. Besides the absolute toler- 
ances and the roughness tolerances, intermediate scale waviness tolerances 
are also specified so that weak shock waves are not produced by locally large 
errors. Existing Langley nozzles are built to very tight absolute tolerances 
(to assure uniform flow) and to very tight roughness tolerances (to delay 
transition). The cost of the 10 inch wide existing Mach 3.5 Pilot nozzle has 
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been estimated to be in the area of $300,000 to $400,000. This nozzle was 
machined from stainless steel and then ground to the shape tolerance. Long 
hours of hand polishing then produced the fined nozzle. Cost reductions had 
to be found to make the quiet- flow nozzles affordable for university research. 
The drawings for the existing Mach 3.5 pilot tunnel nozzle were very kindly 
supplied by Dr. Stephen Wilkinson of NASA Langley, so that design and 
cost comparisons could be made. 

These drawings (see LD-527646) tolerance the area near the nozzle throat 
to what appears to be an absolute accuracy of 0.0003 inches. It seems that 
this tolerance is the reason why the nozzle was ground to shape, at great 
expense. This very close absolute accuracy should result in a nozzle flow 
much more uniform than is usually the case for supersonic wind tunnels. It 
should be remembered that costs increase at least linearly with reductions 
in allowable error; a reduction in this absolute tolerance from 0.001 inches to 
0.0003 inches probably increases costs by more than a factor of 4. It is felt 
that this tolerance is not required for boundary layer instability work, since 
the crucial issue is quiet- flow, not unusually good uniformity. Discussions 
with Ivan Beckwith led to the judgement that this high tolerance might for 
university purposes be applied only to the waviness specification, which does 
relate to flow quietness. Even for the waviness specification it can probably be 
relaxed somewhat. These specifications have been based on Mr. Beckwith’s 
many years of experience. The waviness tolerance is derived from some simple 
computations following from the waviness data found in [4, Figure 4] and [9, 
Figure 4). 

It thus seems reasonable to relax the absolute accuracy requirements to the 
0.001 level of accuracy possible with a quality numerically controlled milling 
machine, thereby reducing the cost of the nozzles by a factor of perhaps 4- 
Waviness tolerances should still be below 0.001 inch per inch, and attention 
will have to be paid to this when the machining strategy is decided on. 
However, the largest error in machining will probably be due to a slow gradual 
variation from end to end caused by misalignment of the milling machine 
ways. This hypothesis was checked by fabricating a test block and measuring 
the machined surface - see Figure 5. This error should have only a limited 
effect on the flow. This change in fabrication technique should by itself bring 
the cost for the nozzle blocks down within budget. 

The material of fabrication is also a crucial issue for machining cost. 
Machining cost for stainless steel is about a factor of 2 more than that for 
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aluminum. However, aluminum does not polish well. Current plans are to 
fabricate the nozzle from easily machined aluminum and then to nickel plate 
the nozzle and polish the nickel. This technique has been used to make x-ray 
optics [5], so there is an extensive literature regarding material selection and 
machining techniques. These optical researchers often use diamond turning 
to produce optics with an axis of symmetry; this may be a cost-effective way 
of maiking an axisymmetric nozzle, since the final surface finish would be 
directly produced without extra finishing. 

As this suggests, the remaining major cost in the fabrication process 
is the cost of polishing the nozzles to reduce roughness. Usual polishing 
specifications are given in terms of the root mean square roughness height 
achieved, which is estimated using various schemes for measuring roughness 
over sample sections of the workpiece. However, for the quiet tunnels the 
crucial issue is the maximum roughness height, which is expected to trip 
the boundary layer locally if it is too large. In a 1986 paper [2] the critical 
roughness Reynolds number (local Reynolds number evaluated at the rough- 
ness height, Rk = puy/p, where all quantities are evaluated at the roughness 
height y = k) was estimated to be between 12 and 42, and the value of 12 was 
chosen for design purposes. The acceptable physical heights depend on flow 
parameters, but generally result in expensive finishes near the limit of those 
normally produced. Normal procedures for machined surfaces involve finish- 
ing using emery cloth, diamond paste, or other abrasives, which produce a 
good average finish as the abrasive size is decreased. However, the quiet flow 
requirement is on the maximum flaw, not on the average finish. Difficulties 
in controlling the maximum flaw have led to specifications on average finish 
nearly ten times tighter than the maximum allowable flaw. Automatic meth- 
ods of measuring the maximum flaw are being sought. Another possibility 
is to find a surface coating which goes on thin and is dominated by surface 
tension while wet. This surface tension would act to smooth out the finish. 
However, it is difficult to find such a coating, since it cannot run during 
application and must sustain reasonable amounts of heat and handling. 

A sample block has been constructed from aluminum using the Purdue 
Central Shop numerical mill. The absolute shape of this block was measured 
at NASA Langley to determine the machining error (see Figure 5). The 
block has since been nickel plated and polished and tested for surface finish 
quality on a computer-controlled stylus machine which also can record video 
microscope images (Tencor alpha-step 200 profilometer with 12 micron sty- 
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lus). Only a small sample of the block could be tested in the machine, which 
is designed for studies of microchip wafers. However, digital records of the 
profilometer traces can be obtained. The finish had a typical peak-peak vari- 
ation of about 7 microinches. The only large flaws were clearly visible, and 
were attributed by the platers to pits introduced during the plating process - 
the largest pit measured had a depth of 100 microinches. Although these pits 
seem less likely to trigger transition thatn equivalent peaks, they must not be 
allowed in the final nozzle finish. The plating shop suggests that it is capable 
of reducing the number of pits drastically if it takes more care (the current 
cost of plating is perhaps 10% of part cost). It would also seem necessary 
to determine a means of inspecting the finished nozzle, and also a rrw^na of 
filling such pits. The good news is that the finishing process in itself does 
not seem to introduce much in the way of large variations - large variations 
from the rms finish will probably be introduced only by the earlier process. 
The sample is still available for tried of various surface finishing methods. 


5 The 2 Inch Purdue Supersonic Wind Tun- 
nel: Modifications for Low Reynolds Num- 
ber Quiet Flow 

During the course of the project development, it became clear that there 
were several issues which could be addressed at very low cost by upgrading 
the current Purdue 2 inch supersonic wind tunnel. This upgrade would allow 
development of instrumentation, calibration of instrumentation, and tests of 
fabrication schemes 5 . The principal requirement was for a vacuum tank to 
allow operation of the tunnel at the very low pressures required to achieve 
quiet flow in an ordinary supersonic wind tunnel test section. However, this 
vacuum tank could then be used for the Ludwieg tube later, as could the 
associated vacuum pump. 

A 500 cubic foot vacuum tank has been procured for this use. A concrete 
paid to support the tank has been installed, and piping to connect the tank to 
the wind tunnel test section has been procured. The tank has been leveled, 

5 Development of high-speed hot film wall sensors is currently being carried out in the 
Purdue 4 inch shock tube, with partial Langley support under NAG-1-1201. These sensors 
are later to be tested in the supersonic wind tunnel, after the vacuum upgrade is complete. 
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and the flapper valve and tee sections fabricated. The Stokes vacuum pump 
has been returned from storage, rebuilt, and reconnected; it is now up and 
running. Unfortunately, there has been an extended delay in hooking up the 
vacuum plumbing, caused by difficulties in the Aeronautics shop (there is at 
present no funding for obtaining outside or Central Shop assistance). We 
nevertheless hope to have the vacuum system up and running by early fall. 

Of course, operation in quiet tunnel mode, even at very low Reynolds 
numbers, will also require settling chamber and nozzle improvements, as well 
as a fine particle filter in the main flow line. A small quiet-flow Ludwieg tube 
seems a more cost-effective method of studying the low Reynolds number 
problem; the improvements described above are of course just as useful for 
such a facility, which has already been proposed. 

6 Summary 

A new kind of short-duration quiet flow wind tunnel has been investigated, 
and preliminary design work indicates that construction should be feasible in 
the university environment. The design of the facility and of the quiet- flow 
test section are described in some detail. 

Although this research area is technically difficult, the author believes a 
suitable program can be carried out in the less sophisticated but also less 
expensive university environment. An inexpensive quiet-tunnel supersonic 
research facility would complement the large and expensive facilities existing 
at NASA Langley, and allow more fundamental research to be done. 
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A 


Large Ludwieg Tube 
Estimate 


Equipment Cost 



Year 1 

Year 2 

Year 3 

L Site Renovation and Building Preparation 




Posts and ties to mount driver tube in lab 

5,000 



Floor rework to enable posts to bear weight of water pressure test 

7,000 



Build working platform and hoist at test section end 

23,000 



Small outbuilding for compressors 


22,000 


subtotal: Site Preparation 

35,000 

22,000 


II. Driver Tube 




48 in. diameter pressure vessel, 200 feet long, with flange for 




contraction end, including welding, pressure test, and all Indiana 




state requirements, for 150 psi working pressure, and for operation 




at up to 225° C } to be erected on site 

60,000 



subtotal: Driver Tube 

60,000 



III. Contraction and Test Section 




Two stage contraction 

25,000 



Pressure containment box 


18,000 


Boundary layer suction system 


7,000 


Supersonic quiet-flow nozzle blocks 

40,000 

75,000 


subtotal: Contraction and Test Section 

65,000 

100,000 


IV. Compressor Air Delivery System 




Air drying system modifications 


3,000 


Oil and particle filters 


2,000 


Pressure relief valves 


1,000 


Pressure gauges, control panel instrumentation, etc 



6,000 

subtotal: Compressed Air System 


6,000 

6,000 

V. Diffuser and Outflow System 




Diffuser sections 


7,000 


Double diaphragm section 

9,000 



Fixed second throat 


4,000 


Flapper valve 


5,000 


Pipe from valves to muffler 



5,000 

Muffler 



5,000 

subtotal: Diffuser and Outflow System 

9,000 

16,000 

10,000 

VI. Model and Instrumentation 



30,000 

Total Direct Costs: 

$169,000 

$144,000 

$46,000 

TOTAL PROJECT DIRECT COSTS (Equipment only): 


359,000 
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B 


Small Ludwieg Tube: Equipment Cost Es- 
timate 


It must be realized that these costs can only be estimated, since this facility 
will be unique. We estimate (Kris Davis, foreman, Purdue central machine 
shop) that a nozzle 14 inches long by 3 inches high by 6 inches wide can 
be built for the following costs. If these estimates prove optimistic during 
preparation of detailed drawings, a narrower test section will certainly be 
feasible. Estimates are as follows: 


1. Cost to fabricate contraction and test section nozzle blocks from alu- 
minum, and to fabricate steel pressure containment box for test section, 
along with flanges. Contours to be cut using numerically controlled 
milling machine. Note that this facility will not have a boundary layer 
bleed slot or a highly polished finish - this is a low Reynolds number 
quiet flow nozzle. Estimate from Purdue Central Machine Shop $15,000 

2. Cost of four 20 foot sections of 12 inch steel pipe for tube: $1400 

3. Pipe fittings for above (elbow and flanges) $1000 

(each weld flange is about $50) 

4. Labor to weld and pressure test the tube and test section $1000 

5. Nuts and bolts for pipe flanges $300 

6. Steel posts to support tube overhead. Estimate includes labor and 

fabrication of connecting hardware. $1000 

7. Cost to fabricate double diaphragm valve section. To be made from 

steel flanges, design adapted from a U. T. Arlington design. Most of 
this cost is Aeronautics shop labor at $25 per hour. $1600 

8. Cost to fabricate wind tunnel diffuser. Current plan is for a square 

design, to be made from steel plate and welded. Includes machining of 
flanges. $2000 

9. Cost to fabricate sliding sleeve section, to allow opening double di- 

aphragm for replacement of diaphragms. Design also follows that of 
U.T. Arlington. $1800 
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10. Cost to fabricate flat plate model for tunnel calibration and receptivity 
experiments. Design to include replaceable leading edge sections. $5000 

11. Total of fabricated equipment estimates: $30,100 
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